PROCEEDINGS 

°F SCIENCE 



O 
O 
(N 

Q 



(N 
> 

O 
(N 



o 



Two-particle Correlation Functions with Distilled 
Propagators 



J. Bulava 

Carnegie Mellon University 



E-mail: |jbulava@cmu . edu 



K. J. Juge' 

University of the Pacific 
E-mail: |k:juge@pacific.edu 



C. J. Morningstar 

Carnegie Mellon University 



E-mail: colin_morningstar@cmu . edu 



M. J. Peardon 

Trinity College Dublin 
E-mail: }n jp@maths . ted . ie 



C. H. Wong 

Carnegie Mellon University 



E-mail: chikhimw@cmu.edu 



For the Hadron Spectrum Collaboration 

Correlation functions of the simplest multi-particle state will be presented using distilled quark 
propagators. The 1 — 2 nn state can be simulated without computing disconnected diagrams and 
thus is the simplest two-particle state that can be studied with quark sources placed on a single 
time-slice alone. We study the quality of the signals of this nn correlation function using the 
quark-smearing guided distillation method. Results will be presented for nn correlation functions 
computed on dynamical, anisotropic lattices. 



The XXVU International Symposium on Lattice Field Theory - LAT2009 
July 26-31 2009 

Peking University, Beijing, China 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



two-particle distillation 



K. J. Juge 



1. Introduction 

Hadron spectroscopy (of excited states) on configurations with ligiit dynamical quarks is reacii- 
ing a stage where one needs to confront the issue of the mixing of single particle states with multi- 
particle states. If the goal is to determine the entire low-lying hadron spectra, one must be able 
to disentagle the single particle states from the multi-particle channels as the thresholds are ap- 
proached with smaller quark masses. 

One of the key issues, then, is the simulation of multi-particle states with controled statistical 
errors. Simulating the excited state spectra without the use of explicit multi-particle operators seem 
to miss some of the excited states due to a poor overlap with the state. One could construct multi- 
particle operators with point-to-all quark propagators, but the need for various momenta operators 
at the source and sink prove to be expensive as is the construction of extended hadron operators. 
All-to-all quark propagators are too expensive for both computational and storage requirements. 
There are various ways of stochastically estimating all-to-all quark propagators in order to over- 
come this problem. In this paper, we test a single-timeslice to all-timeslices propagator method 
which is exact up to an ultraviolet cutoff. 

The simplest two-particle state which can be simulated with this method is the isospin-2 nn 
channel since there are no disconnected diagrams in this channel (see Refs. [p^-[p^] for recent 
dynamical simulations). The method we use is the distillation method (Ref. [|l|]), proposed by the 
Hadron Spectrum Collaboration. The possibility of computing scattering lengths and phase shifts 
via the finite volume method in Euclidean space (Ref. [Q]) will also be discussed as preliminary 
results are presented. 

2. Construction of Operators/Correlators 

The distillation process is described in detail in an earlier paper. We briefly summarize the 
method here for pion and two-pion correlation functions. 

2.1 Distillation 

Hadron correlation functions are usually constructed from quark propagators which have been 
smeared in some way to reduce the high frequency noise from the signal. Jacobi smearing is a 
commonly used smearing method which amounts to a Gaussian-like damping of the high energy 
modes: 



where is the lattice Laplacian operator. Two parameters are tuned to increase the overlap of the 
operator onto the low-energy sector of the particular channel of interest, a and ria. The distillation 
operator on a particular timeslice is constructed from the eigenvectors, Vx'^\t), of the Laplacian on 
that timeslice. 
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where « is a number between 1 and M = NcNxNyN^. We have used n = 64 in this study. We then 
apply this distillation operator on our quark fields which are then used to construct the hadron 
interpolating operators. 

2.1.1 Single Pion Correlation Function 

The correlation function of a single pion operator with arbitrary momentum is given by the 
standard expression, 

In order to simulate pions with definite momenta, we need to sum over the spatial sites at the source 
and the sink with the appropriate phases. 



The sum over the spatial sites makes it necessary to construct the hadron interpolating I/775 y op- 
erator on all spatial points. This is readily done at the sink even with point-to-all propagators, but 
requires substantial more work at the source. The distilled pion operator, on the other hand, re- 
quires no extra work as it contains the propagator from one timeslice to all other timeslices. The 
distilled correlation function is then given by. 
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The eigenvectors and perambulators have been computed and stored for the meson/baryon study 
(Ref. [[n]]). 

2.2 Two Pion Correlation Function 

The two pion correlation function in the isospin-2 channel consists of only two diagrams; one 
of which involves a "quark exchange" ('€') diagram and the other which is simply the product of 
two, single pion correlation functions (Cji{p;t,to)) with back-to-back momenta ('D'). The quark 
exchange diagram is formed by contracting the perambulators in the following way: 
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where repeated indices are summed over. The two pions have back-to-back momenta so that the 
total momentum of the two-particle system is zero. The two pion correlation function is formed by 
taking the difference of the two different contractions to project out the 1 = 2 channel, 

C^n{t,to) = Y,[C4p;t,to)Cn{-p;t,to)-C'j[r{P, -P\tM] 
p 

and the sum over momenta is done to project out the i-wave scattering state. 

3. Simulation/Results 

3.1 Parameters 

We use anisotropic, 2 + 1 dynamical lattices whose tuning of the parameters are discussed in 
Ref. [^. We have results from two different volumes (16^ x 128 and 20^ x 128) with nijj^L 3.5 
and nijiL 4.8 (Ref. [^). The pion mass is roughly 360 MeV on both lattices. The lattice spacing is 
such that ro/us = 3.221(25) with the renormalized anisotropy tuned to = Us/at = 3.5. We use the 
lowest 64 eigenvectors of the Laplacian operator to construct the distilled propagators. The number 
of configurations analyzed were 100 and 94 for the 16^* x 128, 20^ x 128 volumes, respectively. 

We have used pion operators with momenta 0, 1, \/2, and 2 (in spatial lattice units). All 
of the cross correlations were measured to obtain the full matrix of correlation functions. 

3.2 Analysis 

The five-by-five matrix of correlation functions was diagonalized to get the energies of the first 
five nn scattering states. There are several ways of performing the diagonalization, each with its 
own set of advantages and disadvantages. The main concern is the contamination of the signal in 
each of the levels from higher lying states. We provide here some evidence for stability of our data 
against the different methods of diagonalization. The general method of extracting excited state 
energies from a matrix of correlation functions has been laid out in Ref. ||6[ ^. One starts with the 
correlation matrix C,y, 

and solve the generalized eigenvalue problem, 

Cij{t*)w j = X{t* fy)Cij{tQ)w j 

for a given choice of and t* to compute the optimized, fixed-coefficient correlation functions. 
One can then compute the effective mass for each choice and also fit the optimized correlation 
functions to extract the various energy levels. This method does not guarantee that the excited state 
contamination is from states higher than the N'^^ level, but the extraction of the levels does not 
involve any more manipulation of the data and can be fit using a simple, correlated chi-squared fit. 

The excited state contamination can be guaranteed to be from states higher than the N'^ level 
by choosing the ratio of t* to Iq larger than two (Ref. Here we fix this ratio to two and plot the 
effective mass along side the fixed-coefficient method in order to check that the fixed coefficient 
signal has no contamination from the lower-lying excited states. We show examples for the ground 
state and first excited state to show that the systematics are under control (Fig. |l|-^. 

We have also chosen two different sets of (f*,fo) and computed their effective masses to show 
that the dependency there is also very small (Fig. |l|-0). 
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Figure 1: The optimized distilled pi -pi effective 
mass for the ground state (to = 10, t* = 21) on the 

16^ X 128 lattice with 64 eigenvectors. A single 
time-slice was used for the source operator and the 
standard definition of the effective mass was used 
with Af = 1 . 



Figure 2: The optimized distilled pi-pi effective 
mass for the ground state (to = 15, t* = 25) on the 

16^ X 128 lattice with 64 eigenvectors. A single 
time-slice was used for the source operator and the 
standard definition of the effective mass was used 
with Af = 1 . 



3.3 Fitting 

Given the consistency of the various methods of diagonaUzation, we have chosen to fit the 
fixed coefficient, optimized correlators to compute the nn energy levels. One can perform both a 
single exponential fit and two-exponential fits to the optimized correlators, but the two fits agree 
within statistical errors as can be seen in Fig. 5. We therefore extract the two-pion energy from 
single exponential fits whose fitting form was, 



'^-{2M^+SE)t _^^-(2M„+SE)iT-t) 



+ Be 



Here, A,B,Eji;c and M^^ are the constants to be fit and T = 128 is the time extent of the lattice. The 
scattering length is determined from Liischer's formula (Ref. [Q]), 



a, J {atMji){L/as)^ 1 " L/a^ \aj ' {L/usY V«^^ 

Preliminary results for the ratio ao/Mji extracted from uncorrelated fits in units of GeV^^ is shown 
in Fig. 6. The statistical correlation between the single pion correlation functions and two-pion 
correlation function have not been taken into account in the fits, but the statistical error is expected 
to decrease when they are fit simultaneously as they are strongly correlated. Work is underway to 
fit the ratios simultaneously. 
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Figure 3: The optimized distilled pi-pi effective 
mass for the first excited state (to = 10, t* = 21) on 
the 16^ X 128 lattice with 64 eigenvectors. A single 
time-slice was used for the source operator and the 
standard definition of the effective mass was used 
with Ar = 1 . 



Figure 4: The optimized distilled pi-pi effective 
mass for the first excited state (to = 15, t* — 25) on 

the 16^ X 128 lattice with 64 eigenvectors. A single 
time-slice was used for the source operator and the 
standard definition of the effective mass was used 
with Af = 1 . 



4. Summary 

The simulation of the I = 2 nn scattering states using distilled quark propagators on 2 + 
1 anisotropic, dynamical lattices have been presented. A relatively light quark mass was used 
on two different volumes to check for finite size effects. The distilled propagators have allowed 
the extraction of the 4'^ excited state to 10% errors on 100 configurations. Work is underway 
to incorporate the correlations between the single pion and the two-pion measurements to reduce 
the statistical errors in the determination of the scattering length and to determine the scattering 
phase shift. A modified version of the distillation procedure is also being investigated to reduce 
the computational/storage cost of the method in larger volumes and correlation functions which 
contain disconnected diagrams. 
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Figure 6: Preliminary values for the scattering 
length in units of GeV^^ extracted using the 1/L 

Figure 5: The f,„,„ plot of the fitted values for the ^^^^^^ ^^he l/L^ value is shown for ref- 

pi-pi energies of the second excited state (V = 16 ). 

The dashed fit values indicate a poor fit. 
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